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Abstract 

Co-sputtered Au-Ni thin films having thickness of 30 nm were deposited on Si(100) 
substrates and irradiated with 160 keV 40 Ar + under ambient condition at a number 
of fluences and analyzed using Rutherford backscattering spectrometry (RBS). The 
variation of Au signal counts in the RBS spectra with ion dose has been investigated. 
The distribution of Au, Ni and Si atoms over various depths within the as deposited 
and irradiated samples have been computed using the backscattering data by means 
of a direct analytical method. Au and Si profiles have been fitted with error function 
and the relative changes in variance for various ion fluences compared to that of as 
deposited profiles have been studied. The spreading rates of different constituents 
across the interface due to Ar ion impact have also been discussed. 
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1 Introduction 



As the dimensions of solid state devices are scaled down, it is becoming in- 
creasingly necessary to improve the performance of the interconnects and their 
electrical characteristics. Metal silicides have some advantages over pure metal 
contacts due to their low-electrical resistivity and good thermal stability and 
these silicides are extensively used in microelectronic devices as contact mate- 
rials to the source/drain areas to control the Schottky barrier height [1-3]. En- 
ergetic ion beams, when penetrate through the interface of different materials, 
produce massive atomic transport across the interface which results in many 
stable, unstable or even thermodynamically non-equilibrium phase formation 
around the interface. Due to improved electrical and chemical properties of 
the ion irradiated materials, mixing of metal-metal or metal-semiconductor 
systems using ion beam is frequently used to tailor different suitable contact 
materials for electronic devices [4,5]. Ion bombardment induced composite 
formation in Au-Ni bilayer or multilayer films deposited on Si [6, 7] or III- 
V nitride semiconductors [8] have been proved to be an effective method to 
fabricate contacts having low resistivity (typically, ~ 6 x 1CT 4 fi-cm [9]). 

Rutherford backscattering (RBS) technique is an easy and efficient tool to 
study the dynamics and atomic rearrangement around the interface. Analy- 
sis of RBS data and extraction of elemental depth profiles are mostly done 
using some standard codes [10-12] which are based on matching of experi- 
mental spectra with the theoretically simulated spectra of hypothetical sam- 
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pie structures having multiple layers of different compositions. This type of 
indirect depth profiling is a difficult and time-consuming job sometimes leads 
to ambiguous results. Even for RBS spectra showing well separated peaks 
of constituent elements, this type of indirect method to synthesize atomic 
depth profiles is not always a trivial issue specially when the number of layers 
in the hypothetical sample is large or when one investigates the accumula- 
tion/depletion of a particular element over a small region inside the sample. 

The aim of the present study is to implement a direct method to extract the 
depth profiles of the constituents from the Rutherford backscattering spectra 
and to get insight of the interactions between thin co-sputtered Au-Ni film and 
Si substrate during ion bombardment. The extracted concentration profiles of 
such a system and its analysis allows investigation of the mass transport of ma- 
terial across the interface during Ar ion irradiation, which may be important 
in the field of interconnects and contacts in the electronic device fabrication. 
Specially, relatively lower concentration of Ni in the system introduces control 
over the interface dynamics and formation of Au(Ni)Si alloy generated due to 
ion irradiation. In this study, we have shown the nature of impurity spreading 
into the matrix and subsequent broadening of interface due to 160 keV Ar + 
ion impact. 

2 Experimental details 

Thin metallic films of co-sputtered Au and Ni were grown on cleaned and 
polished Si(100) substrates in Ar atmosphere by dc magnetron sputtering. 
Prior to the deposition, the Si substrates were etched using HF solution to 
remove the native oxide layer from the surface. Before Ar introduction in 



3 



the deposition chamber to create plasma necessary for sputtering the base 
pressure of the chamber was 8 x 1CT 7 mbar. The deposition was carried out 
at room temperature and the pressure of the chamber during deposition was 
9 x 1CT 3 mbar. During deposition, both the Au and Ni sources were operated 
simultaneously with same power to get uniformly mixed single layer of Au and 
Ni over the substrates. However, different sputtering rates of these elements 
gave rise to unequal stoichiometric composition within the films with 90 at% 
Au and 10 at % Ni. Thickness of the deposited layer was ~ 30 nm as measured 
using x-ray reflectivity technique. The deposited samples were then irradiated 
with 160 keV 40 Ar + ions using high current ion implanter (Danfysik) facility 
of Saha Institute of Nuclear Physics. A schematic diagram and details of the 
facility can be found elsewhere [13]. The fluences of the implanted ions were 
varied from 5 x 10 14 to 5 x 10 16 ions/cm 2 and the pressure in the target 
chamber was 3 x 10~ 7 mbar. Homogeneous irradiation was achieved by means 
of a magnetic X- Y beam scanning system. No special arrangement was made 
to control the temperature of the samples during implantation but the beam 
current was kept low to avoid beam heating induced atomic transport within 
the sample. Depending on the attainable ion fluences, the ion beam fluxes were 
kept around 400 nA/cm 2 . The projected range of Ar ions in the target was 
determined using SRIM2006 code. The irradiated samples were analyzed by 
means of Rutherford backscattering spectrometry (RBS) using 2.05 MeV He 2+ 
ions obtained from the Pelletron of Institute of Physics, Bhubaneswar. The 
backscattered data was collected using a silicon surface barrier detector with 
16 keV resolution placed at 160°. Corresponding depth profiles were extracted 
from the RBS spectra by direct method described in the following section. 
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3 Scheme of data analysis 

The depth profiles of the constituent elements have been extracted directly 
from the backscattering data by a computer program written using MATLAB 
package. The method is particularly suitable for this system as the peaks 
corresponding to different elements are separated and well resolved. The di- 
rect determination of stoichiometry of the ion irradiation induced atomically 
mixed layer ensures accurate determination of elemental composition at vari- 
ous depths compared to other indirect simulation based standard codes. How- 
ever, using both the above-mentioned procedures for a given analysis could be 
very convenient and helpful to avoid misinterpretation of the results. 

Under this data analysis scheme, different elemental edges for the incident 
He 2+ ion energy E Q , which is 2.05 MeV in this case, are calculated and tracked 
in the experimental spectra using the relation k x Eo, where k x s are the kine- 
matic factors with x representing the major elements present in the sample, 
namely, Au, Ni and Si. The signal height (H x (k x E )) of an element x (= Au, 
Ni or Si) at its respective edge is given by [15] 

H x (k x E ) = a x (k x E )QnN™ x T x /cos9 (1) 

where a x (k x E ) is the differential Rutherford scattering cross sections of el- 
ement x at energy k x E , Q is total amount of charge accumulated on the 
sample due to He 2+ incidence during the backscattering experiment, fl is the 
solid angle subtends by the solid-state detector at the sample and 9 is the 
angle subtends by the incident He beam with the surface normal (in the cur- 
rent setup 8 = 0°). N™ lx is the atomic density of element x in the mixture of 
Au, Ni and Si at the surface and t x is the thickness of the imaginary surface 
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layer inside the sample which is so chosen that the energies of He particles 
backscattered from this layer after colliding with atoms of element x fall within 
the energy window k x E — £, where £ is the energy width of one channel of 
the detector. It is evident from the above fact that the depths r for different 
constituents will be different depending on the concentration of a particu- 
lar element in that mixture, respective kinematic factor and scattering cross 
section. The atomic density of the mixture N mix is given by the relation 

N mix = mN Au + nN Nl + pN Sl (2) 



where iV's are the atomic densities of elements indicated by the subscripts and 
m, n and p are relative fractions of Au, Ni and Si respectively. If the presence 
of the implanted ions is neglected, the sum of all the atomic fraction is unity 
(i.e. m+n+p = 1) and in terms of relative fraction and atomic density, N™ lx 
can be expressed as 

NZ x = mN Au 

N™ x =nN m (3) 
N£ ix = pN Si 



Using equation ([3j), the signal heights of the different elements present on the 
surface (equation fll}) can be modified as [15] 

H Au (k Au Eo) = (J Au (k Au E )QVLmN Au T Au (4) 
H Ni (k Ni E ) = (T Ni (k Ni E )QQnN Ni T Ni (5) 
H Si (k Si E ) = crsi(ksiE )QflpNsiT S i (6) 



The thicknesses corresponding to different elements expressed by t x are related 
to the stopping cross section factors of respective elements ([e(Eq)]™ ix ) by the 
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expression [14] 

t x = — i — (7) 



and the weighted sum of r's gives the effective thickness of the layer: 

r = mr Au + nr Ni + pr Si (8) 

Therefore, equations (HJ)- (JBJ) can be rewritten in terms of thickness and stop- 
ping cross section factor as 

H Au (k Au E ) = a Au (k Au E )QVtmN Au ^ . (9) 

[e(E )] A l N mix 

and so on. 

Using the above mentioned expressions of signal heights of elements present 
within the surface layer r, the ratio of relative fractions can be expressed as 

m _ H Au (k Au E Q ) a Ni (E ) N m {e(E Q )] A u ^ 



n H Ni (k Ni E ) a Au (E ) N Au [s(E )] 



mix 
Ni 



(10) 



and 



m = H Au (k Au E ) a St (E ) N St [e(£ )g 
p H S i(k St Eo) a Au (E ) N Au [e(E )]™ x 
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The ratios [e(^C/[«]r and H E o)\Z I 

™ x can be treated as 

unity in the first approximation and using the tabulated values of cr's [15], 
relative concentration of each element has been estimated. 

As the concentrations of all the elements within the depth r can be calcu- 
lated using the above approximation, the values of [e^Eq)]™* / [e(Eq)]™* and 



7 



[e(-Eo)]™" / [ £ (E )]™. X can be extracted using the expression of stopping cross 
section factor of element x: 



he,, 



-\mix 



k 1 

n x mix I 771 \ I x mix 1 ;„ 771 
COS COS ffl 
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where (= 20°, in this experimental setup) is the complement of scattering 
angle and e mix (Eo) and e mix (k x Eo) are the stopping cross sections of the mix- 
ture at energies E and k x E respectively. The expressions of stopping cross 
section of the mixture for any given energy E is represented by 

e mx (E)=me Au (E)+ne m (E)+pe S i(E) (13) 

where e Au (E), e m {E) and esi(E) are individual stopping cross sections of Au, 
Ni and Si respectively which are evaluated by the polynomial fit expressions 
[15] 

s(E) = j2A i E i (14) 

i=0 

with A — A 5 being constants for a particular element when E is expressed in 
MeV. 

Using equations (fT2l) and (jTBl) . the stopping cross section factors for different 
elements are calculated and fed in equations (TTUT) and ( Till to calculate ex- 
act values of atomic fractions which values are again used in equations ffl2l 
and (fTBl in iterative manner to produce more accurate values of elemental 
concentrations in the layer. 

As the thickness (r) and composition ^N™ 1 ™) of the first layer is known, the 
energy, loses by the penetrating He ion while traveling through that layer can 
also be estimated and the energy of the ion just before entering the next layer 
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(E) is given by 

e = e' — — [f^jjy — 

COS# 



where £" represents the energy of the ion at the point of entering into the 
previous layer (= E , when dealing with second layer starting from surface). 
Simultaneously, signal heights of each element in the next lower channels start- 
ing from the corresponding elemental edges are also extracted. The ratios of 
atomic fractions in the second layer are given by 



m = H Au (E Au ) a m (E) N Ni f [e(E)]™ e(E m ) e(k Nl E) 
n H m (E m ) a Au (E) N Au \ [e(E)]^f e(E Au ) s(k Au E) 



(16) 



and 



m _ H Au {E Au ) a Sl {E) N Sl f [e(E)]™ e(E Si ) e{k Si E) 



(17) 



p H Sl (E Sl ) cr Au (E) N Au { [e(E)}™ x s(E Au ) e(k Au E) J 

where H Au (E Au ) is the height of Au signal in the channel preceding to that 
which corresponds to Au edge and E Au is the energy equivalent of that par- 
ticular channel. Similarly, Hjsn{E^j) and Hsi{Esi) is the heights for Ni and Si 
signals respectively at channels next to Ni and Si edges towards the lower en- 
ergy side. e(E Au ) and e{k Au E) etc. are stopping cross sections at energies E Au 
and k Au E Au and so on. [e(E)]™ tx is the stopping cross section factor for ele- 
ment x present in the mixture and expressed by the same relation as equation 
f fT2l) with energy E replaced by E. 

When treating the second layer, the approximation applied to the surface 
layer for determining the elemental composition is also valid which treats 
the values of the quantities within curly braces in equations (jTBI) and (flTI) 
as unity. So, first estimation of the concentrations can be figured out using 
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that approximation which, in turn, can be used along with equation (fT4|) and 
tabulated values of a's [15] in iterative manner to determine accurate values of 
elemental concentrations in that layer. Subsequent determination of N mix and 
r for the layer can be done using equations (T5]), ([7]) and ([8]). For estimating 
elemental concentrations of third, forth and subsequent deeper layers, same 
formalism and expressions as applied to second layer can be followed. 

It should be noted that experimentally obtained RBS spectra are convoluted 
with the apparatus function which contains initial beam energy dispersion, 
energy resolution of the detector, energy loss and energy-loss straggling in 
the detector as well as in the sample. But the determination of deconvoluted 
spectra using the normal deconvolution process is not trivial and sometimes 
introduces unwanted undulations in the spectra which may yield negative con- 
centration of an element at some depth. So some other approaches described 
in ref. [16] or Bayesian probability theory [17] are suitable for recovering the 
'true' spectra. Since the final aim of this paper is to investigate the relative 
broadening of Au and Si profiles with increased ion fluence in case of co- 
sputtered Au-Ni system, the effect of apparatus function causes unidirectional 
change in the RBS spectra and hence nullified. 

4 Results and discussion 

Figure 1 represents Rutherford backscattering spectra of as deposited and 
samples irradiated with Ar ions at different fluences. It can be observed from 
the figure that the Au signal height decreases with increasing dose which 
signifies sputter erosion of the films. On the other hand, Ni signals do not 
decrease significantly with increasing dose. This is due to the fact that higher 
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sputtering yield of gold triggers preferential sputtering of Au atoms within the 
mixture of Au and Ni. As a result, Ni signal heights decrease less rapidly than 
that of Au. Along with diminution of Au and Ni peaks, tailing of these signals 
towards lower energy side can also be observed representing ion irradiation 
induced atomic diffusion of constitutes in the direction of sample interior. 
The accumulation of Ar atoms which penetrate into the samples during ion 
irradiation, show no considerable change in the spectra until a dose of 1 x 
10 16 ions/cm 2 beyond which a strong peak of Ar emerged between Ni and Si 
peaks. In contrast to the spreading of Au and Ni signals towards Si side, Ar 
peak at a fluence of 5 x 10 16 ions/cm 2 shows symmetry around its maximum 
position implying Gaussian type distribution commonly observed in case of 
implanted ions in solids at relatively low doses. 

Variation of total counts enveloped by the Au signals in the RBS spectra with 
ion dose have been showed in figure 2. The total count for as deposited sample 
has been taken as unity and other counts are normalized with respect to that. 
The counts decrease rapidly at the initial stage of ion irradiation and after 
achieving a dose of 1 x 10 16 ions/cm 2 it becomes monotonous which is in 
accordance with the previously reported result regarding sputter erosion of 
thin films [18]. The data has been fitted using exponential decay function 



where C and Cq are the normalized and saturated counts respectively and 
A, B are constants. The fit gives the values of different parameters as C = 
0.62 ± 0.03, A = 0.36 ± 0.04 and B = (27.98 ± 11.09) x 10 14 ions/cm 2 . The 
exponential decay of the Au counts with ion dose suggests that sputtering of 
thin films is different from that of bulk material. At the early phase of ion 




(18) 
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bombardment the sputtering is to some extent analogous to bulk sputtering 
and the sputtering yields of Au and Ni for an ion-target combination as the 
present one are ~ 4.0 and ~ 0.5 atoms/ion respectively, as calculated using 
SRIM code. The outcome of larger yield for Au is a Ni-rich top layer along 
with surfacing of Si due to ion induced cascade mixing of the components 
which produces a saturated region in figure 20 above a critical dose of 1 x 
10 16 ions/cm 2 . 

The depth profiles of as deposited and irradiated samples as shown in figure 
3 reveal the possible cause of exponential type decaying of integrated Au 
counts represented in figure 2. The depth profiles of as deposited and sample 
irradiated at a dose of 5 x 10 14 ions/cm 2 show no significant difference with a 
presence of ~ 90% Au and ~ 10% Ni at the surface. After achieving a fluence 
of 1 x 10 15 ions/cm 2 , the sputter erosion of top surface layer is visible in the 
Au profile which is more noticeable at the next dose by the advancement of 
Si profile towards the surface. The Si concentration at the surface has reached 
around 15% for ion dose of 5 x 10 15 ions/cm 2 , with a reduction of the amount 
of Au. At higher doses the effect of preceding of Si profile is more dominant. 
Moreover, at a fluence of 1 x 10 16 ions/cm 2 , a surface layer of 10 nm having an 
approximate composition of Au5sNi 5 Si37 covering another layer of Au4sNi 5 Si5o 
which spreads over ~ 20 nm has emerged. At the next higher ion dose, i.e. 
at 5 x 10 16 ions/cm 2 , the buried compound layer has been exposed leaving 
completely mixed surface layer having a thickness of ~ 20 nm. The formation 
of mixed layers at the surface and subsurface regions at higher irradiation 
doses are responsible for saturation of Au signal count represented in figure 2. 
However, the Ni profiles at various fluences show very little observable changes 
due to irradiation. 
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The extracted atomic profiles has been used to quantify the growth of the 
mixed region due to ion irradiation. Au and Si profiles are fitted with error 
function to deduce the change in variance a 2 as a function of ion dose. The 
relative change in variances (Act 2 ) for various doses with respect to the vari- 
ance of Au profile for as deposited sample is represented in figure 4. A linear 
growth of the variance with the ion fluence had been observed previously in 
many metal/semiconductor systems due to ion irradiation [19-22]. But for the 
present system, the linear variation of Ac 2 with ion dose seems to be valid 
only up to a fluence of 1 x 10 16 ions/cm 2 . Beyond that Acr 2 for both Au and 
Si profiles decrease to a lower value. The linear fit of present data up to 1 
x 10 16 ions/cm 2 give the mixing rates Acr 2 /fluence = 6.88 ± 0.77 nm 4 for 
Au and 9.34 ± 0.43 nm 4 for Si. The difference in mixing rates is probably 
due to the diffusion of Ni atoms in Si layer which occurs simultaneously with 
Au. The lower heat of compound formation (-16 kJ/g at.) for Ni and Si and 
hence ready creation of different Ni-Si phases forms a physical barrier before 
the Au atoms to spread at an equal rate of Si. The decrease in Aa 2 beyond 1 
x 10 16 ions/cm 2 can be explained with the help of figure 2 and 3. The SRIM 
simulated range of impinging Ar ions is ~ 125 ± 65 nm which is sufficient 
to deposit considerable energy on both sides of the interface and produces 
collision cascades. Collision cascades both in film and substrate regions set 
off recoil implantation, cascade mixing and radiation enhanced diffusion pro- 
cesses. Simultaneous preferential surface sputtering and cascade mixing not 
only form a Au-depleted surface but also an accumulation of Si around the 
near-surface region at and beyond 1 x 10 16 ions/cm 2 Ar dose (figures 3(e) and 
3(f)). Eventually, the evolution of homogeneously mixed layer of Au45Ni 5 Si 5 o 
which has been already discussed, contributes towards the sputtering process. 
The notable fact about these mixed layers is the absence of preferential sput- 
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tering of any species due to ion impact. The major components of the mixture, 
namely Au and Si are sputtered at the same rates which is supported by the 
SRIM simulation of similar kind of target composition-ion combination. The 
simulation provides the erosion yields of both Au and Si as ~ 1.5 atoms/ion 
and as a result, not any single component but the whole surface is sputtered 
away with ion impact. This fact is further supported by figure 2 which shows 
a flat region past a threshold dose. Thus, the bulk Si is becoming exposed 
which gives rise to lower values of Aa 2 at an ion fluence of 5 x 10 16 ions/cm 2 . 

5 Conclusion 

In this study, an efficient analytical method to analyze RBS data has been 
presented. Using that method, depth profiles of different components of co- 
sputtered and ion irradiated Au-Ni thin films on Si has been investigated. It 
has been observed that variation of Aa 2 with dose follows a linear relationship 
up to a certain ion fluence which yields mixing rates of 6.88 ± 0.77 and 9.34 
± 0.43 nm 4 for Au and Si, respectively. The lowering of Aa 2 values ahead 
of a critical dose can be explained by the formation of Au-Ni-Si composite 
at the surface from which elements are sputtered out evenly which is further 
supported by variation of Au counts in RBS spectra with ion dose. 
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Figure captions 

Figure 1: RBS spectra of (a) as deposited co-sputtered Au-Ni sample and 
samples irradiated using 160 keV 40 Ar + with doses of (b) 5 x 10 14 , (c) 1 x 
10 15 , (d) 5 x 10 15 , (e) 1 x 10 16 and (f) 5 x 10 16 ions/cm 2 . The upper and 
lower panels represent the same spectra but in the lower panel each spectrum 
is vertically shifted from others. The spectra are normalized with respect to Si 
signal height and different scales are used for Au and Ni/Si signals for clarity. 

Figure 2: Integrated counts under Au signals in the RBS spectra of co- 
sputtered films. The counts are normalized with respect to that of as deposited 
sample. The data has been fitted using the equation of exponential decay. 

Figure 3: Depth profiles of (a) as deposited co-sputtered Au-Ni sample and 
samples irradiated using 160 keV 40 Ar + with doses of (b) 5 x 10 14 , (c) 1 x 
10 15 , (d) 5 x 10 15 , (e) 1 x 10 16 and (f) 5 x 10 16 ions/cm 2 . 

Figure 4: Change of variance of (a) Au and (b) Si profiles as a function of 
ion dose and linear fit (dotted line) of the experimental data. 
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